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ABSTRACT 

The purpose of this investigution it to uetermine the 
effect on the temperatures anu the perforgunce peraneters of 
tne addition of (1), liquid hydrogen, anu (2), liguid ammonia 
to three bipropelieant systems currently of interest, These 
bipropeliant systems sre: nitrogen tetroxiae-hycrazine, 
hydrogen perautée-ayéresins, atu RPNA-hycrazine, 

Por each of the six tripropellant systems investigated 
the chamber temperature, exhaust teaperature, and ull tiie 
important parameters were calculatec, These calculstions 
for each system were carriec out for both equilibrium flow 
anu constant cosposition flow assumptions, All the results 
are listeu in tables @nu the Bore imsportunt paraweters ere 
presented in graphicai form as well, 

One hyurogen containing system am one ssmonia con- 
taining system are ciscusesed in detail, General remuite 
for all systems ineluce the fact thet the largest proportion- 
ate aecreases are observec in the varixtion of the exhaust 
ana chamber temiperutures wiereas the specific ifpulse shows 
slight increasws or proportionate uecreuses of much saller 
magnitude than the proportionate decresses in temperstures, 

veveral New perameters are introducei in an eftort 
to Bake possible the prediction of performances of « tripro- 
pellant system ir the chamber temperature variation is known, 
These pardweters may s1S80 be usec in un inaireet spanner to 
illustrate the relutive merit of liquiu hydrogen anc liquid 


ammonia as coolants in tripropelliunt systexs, 
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INTRODUCTION 


One ol the ever present problems in the design of 
any heat engine is the problem of removing the heat generated 
by the combustion at a rate hizh enough to prevent failure 
of the engine, In a thermal jet or rochet motor the temperae 
tures of combustion and rates of heat transfer are greater 
than in any other practic:l heat engine, It is therefore 
easy to understand why this problem is one of the major 
problems in thermal jet and rocket motor design, 

One solution to the problem is regenerative cooling. 
A regenerative cooling system maintalius the motor wall ma- 
terial below its criticul temperature by the circulation 
or the fuel or oxidant outside the motor wall, For many 
of the high performing propellant systems this method is 
inadequate and more efrective cooling methods must be used 
(cf, Ref, 1), Furthermore, many of the newer fuels and 
oxidants are not stable enouzh to permit their use in rea 
generative cooling systems, 

Theoreticaily at all operating teuwperatures, and 
practically, at lower temperetures, it hes been shown that 
film or transpiration cooling systems have certain advantages 
(Gf, Ref, 2). Also it appears that theoretically such cool- 
ing techniques are applicable to very hot propellant systems 
over longer periods of tine, 

In film or transpiration cooling, the coolant moves 


against the flow of heat through the wall material. On the 


hot side of the wall the coolant elther evapor:tes, or 
tanes part in the chumber reaction, or some combination 
of both, Thre extent of the mixing with the combustion 
gases is not yet known nor is it known what efficiency 
ot absorption of rediant energy by & gas or Vapor streas 
may be expected, | 

If adeyguete mixing eno complete chemical] rection 
are aSsuped, the aduition es « film or trenspiration coolent 
of a properly chosen third component Bay have certein de- 
sirable effects on the calculsted performenc® 4m edinabatic 
flume tewperature, aAlthoug» tiw extent to enich che os 
Sumptions are velidg must be Getwurmine es, arleentelly, 
the theoretical effects on peri ore. oce anu gas lebper ture 
are of impediate interest since they will indteste the at- 
rection anu orcer ot magnitude of the effects thet can be 
expectwa (Cf, Nef, 9), 

Although liquid water has bewn seriowsly con» ered 
(Cf. Refs, 7 anc S), it has bewn previousl, Shown that the 
most efrect.ve thirc components are liquia or gaseous hy- 
drogen aud liquid or gaseous ammonia (Cf, Aefs, 5, 4, 6 ana 
7). These coolants or their cissouistion products lower 
the chamber temperiture by aecreesing the available energy 
per pouna of total propellant, The performance may increase 
because of the lower average molecular weight of the pro- 
duets of reaction. If, however, the performance cecrenses 
its proportionsete decrease is oslways such less than the 


proportionate decrease in chamber temperature, 


Fit, 
wee) ed 


For bimes ence the investigation af six tri- 
propéllant systems currently of interest was undertaken 
and the results are presented herein, The systems investi- 
gated ares nitrogen tetroxide-hydrazine with liquid hydrogen 
wedded, nitrogen tetrcoxide-hydrazine with liquid ammonia 
added, hydrogen peroxide-hydrazine with liquic hydrogen 
added, hydrogen peroxide-hydrazine with liquid ammonia 
added, RFNA-hydrazine with liquid hydrogen added, and 
RFNA—hydrazine with liquid ammonia added, 

fhe performance parameters were evaluated for each 
triprovellant system at a mixture ratio corresponding to 
stoichionetric proportions with respect to the bipropellant 
oxidant and fuel, but with varying amounts of liquic hydro- 
gen or liquid ammonia added in excess. Stoichiometric pro- 
portions of the bipropellunt components were chosen since 
the chamber temperature is very nearly the maximum value 
at this mixture ratio and the reduction in flame tempera- 


ture resulting from an added component would be most marked, 
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thuber of eeiys oF wutwe waper (10) tn 
che prouucta of reaction 
products of reaction 
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Musuor of moles of atseic hycregen ome 
an the produets of rvuction a 

fuaber of solus of caygut (uy) in ule 


products of reection - 


Wunber oF moles oF atowic oxyutn (G) im 
the prowucts of re.etion 

Wuaber of soles of nitrous oxide (NO) in 
tw prouucts of resection : 

Wusber of poles Of nitrogen (My) tt the 
proubets of rection 

| Wuliber of gram Atoms of hydrogen in the 
react nts 

Rusder of gheb atoms OF nitrogen in the 
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reactsnts 

Wuwber of groom atoms of Oxygen in the 
Pewctalits 

Velocity of sew corresporming te tine- 
bur cOnuitlots (ft see~+) 


Ereective €2b.0st veLoeity (ft vac-t) 


C*¥ 


Cr 


(Cp), 


(Cy) 


C2 


“(Isp) 


L(y 
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ht... 


—_— . ~1 
Charéeberistie velocity (@t sec ) 


Taéoretieali thrust ceefrlesent of nezete 
Apparent molur isobdaPle Reat capacity of 
products of rézvetion at @fudii Prius 
chamber temperature (cal not ye) 

Apparent molnr isochoric heat Cepacit; 
o© products of re&z@tion .t Suit ipriam 
chamber tempernxture (cal moi7+ ox-t) 

Averayve apparent molar isobsric heat 
Cipecity of products of reaction during 
their passage throuxgh nozzle 


1 o,-4) 


‘(Gud Goal - 

AVerage upparent Bolsr isochoric heat 
capacity of preauucts of reaction cuPing 
bivilY passage through nozzle 


(cal molt O,-1) 


Dimetislonless parameter equal to 


if « & 
isp = Te 


Dimensionless puremetéi Eiual to 


a -% 
h? te” 


Nozzle throat area (sq in) 

Thrust of rocnet motor (1b) 

Acceleration of gravity (arbitrarily 
chosen eyual to 3.,2 ft seen”) 


Altituae index (nmi) 





he 


NHS), am’, etc, Clumge in wutiaipg of the prowucw of 
rtection (» cil) Detwotn superseript 
veryeratare ([%) = stboerlyt ten 
per ture (72) 

ntheiys clange ip » spwe. fied aysten 
froe T, to Ty ith the syatew in etepa- 
ewl @yublioriugw (& c#l) 

Specific impulse of propellant (sec) 

Mechific.i @juivesent of heat (&,lee 
x 1% ergs & cain) 

By by, etc, Byulliorius constents for pert.cular 

; reuctious ¢apressw in terns of partial 

pressures of components - 

Any» Any» etc, Squilibriue conStents for particulier 

> Presctious #apresce. in terea of nueber 
of moles of componcats 

TOtel] weight of recctunts (gH) 

Weight rate of flow through nozsi« 

(ib sec™*) 

Average molecular ee@igit of the provwucts 
of reaction uring thelr paseese through 
tre nos.le 

Average Boleculer welgi.t Of the prowucts 
of réection at eyuilldriusy clasver 
tegpere ture 


Numoer of moles of predutts of resection 
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oe etc, 


Am, Atm, FC. 
Fue 


P 
Pe 
Po 
Te 2000 
Yay? ey ;, ete. 


fats 
av Te 


dp (reactants) 
Ip (preducts) 


RNA 


» 
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AVcraxe number of moles oF prouucts of 

reaction pre.lent during their passage 
throuzn the nozile 

Nunber of Holes at temperature (°K) 
indicated by superserigt 

Humber of moles of proaucts minus the 
number of wolcs cf reactants inaicated 
in baiancea equiisibriun equation ror 
a particular reuction 

Total pressure of the prouucts of 
reaction (atn,) 


Noz.lé exhaust prevoure (psia) 


Chai.ber pressure (psia) 


Heat aviilab.ie Crom completion of re- 
action at temper.ture (°K) as indi- 
cateu by superseript (k cal) 

Change in heat availebie between T, 
Athi Ls Cn eal). 

Heat of formition of the renctants @t 
zgo? kh (Ck eal) 

Heat of Yform.tion of the prouuets at 
300° kK (k cz#1) 

Red fuming nitric acid (in this investi- 
gation, nitric acid with 6.6% by 
weight NO.) 

Universal gus constant (1.666 cai 


nio17 1 7, 6810 & 10° eres wont oym) 





7, Bra? 








= 5 


Absolute tecperature (C8) 

Equilibrius chaeber tewperature (°%) 

Paheust tewpereture (°F) 

Averugt density of rocket propellant 
(gm cn-*) 

Ratio ef isovaric to agecnoric hust 
Capecity (C,/Cy) 

Ratio of apperent isobaric to isechoric 
heat capacity of products of resection 
@t eguiiibrius clhesber tecpersture — 

Ratio of average Bpparent isobaric to 
isochoric heat eapacity of the pro- 
ducts of reectiou during their paasage 
through tle nozile 

A function of § defimea by 
rhe yep ty 


ri3 evaluated for F = 


_ 


c 
Superscript zero inulc:tes value of 


parameter for bl,propelleut system 


ot SboiGuicketric sLxtvure ratio 


The results et this investigation can be better 
evaluatea when they are comsidereu from the point of vies ef 


=> 
tne assumptions on which they are based. Hence, bercre the 


methods of 
lustrated, 
concerning 
prevail ir 
ly valid, 


fhe assuuptions are listed below in the approximate 
order in which they ere appliea in Parts Il ana Til, k-ny 


of tise assumptions are uiseusseu in greater uet il .in these 


suceeéding 


i. 


Cunleuluting the results are uiseussec anu ile 
it is necessary to set tiown the assumptions mace 
the combustion ani Plow conditions which wust 


the theory ana @eyguations used are to be comnplete- 


sections, 


OF ADSUAPTIONG ARU IRLTRMUDUCTION OF PAKABETARS 
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First of all, it is assumed that the pro,;ellants. 
reacl completely ana there is sufficient time , 
for establishuent of eguillbrium concentrations 
of Une normally unméxcite. cornponents (Ho, Hod, 
OH, BH, Ogg 0, NO, No, W@W) at th adiwbatic flume 
temperature determinea by the mass and neat 
balance, 

All iwsinor co#ipenents sre considere. in these 
calculations except atosic’ nitrogen (N), Neg- 
Lecting this cotponment introduces an error 
which is smuller thun possible errors in the 
equilibrium constants which uvetermine the 


composition, 


Ss 


~10- 


The equipertition of energy among the electronic, 
vVibrationzl, anu rotstions] energy levels is instan- 
taneous, vot. in the chamber and in tne nozzle ex- 
pansion processes for thew constent compositien and 
equilibrium flow comitions assumed, 

The gus Beicture composition, dO Hib, 3 nn ate 

are assusec to be linesr over température ranges 

of one hunureu degrees Kelvin, 

For calculution of the characteristic velocity, the 
Fatio of specific heats (y¥) is assumed to be tie 
ratio of the specific heats at the equilibrius chas- 
ber temperature (%.). This was caleuleted on the 
basis of the enthelpy change over the one hundre« 
degree temperature interve] nearest the chamber tem- 
perature at constant pressure, the chamber pressure, 
Por calculation of the exhaust temperature, a con- 
stant average ratio of specific heats (y) is used, 
This was calculated on the basis of the enthalpy 
change during the mozzle expansion process, 

The function !, defined below, can be represented 
by a linear function of ¥ (Cf, Khef, 5), By 
definition, 


[F 26 ay FO 


By linear approximation, 
r= ~1047 # ,S068¥ 
This linear approximation was checked over the r rge 


of ¥ encountered in these celculatiorns ana found 


10, 


li. 


lc. 


i 


to be @qual to the vilue cxiculatea by the Longer 
expression to within one-half per cent, 

In the calculations in which equilibrius flow is 
assumed to be maintained in the nozzle expansion pro- 


cess, the gas mixture composition changes continuous- 


ly and instantaneously in accordance with the tempera- 


ture and pressure variations as determined by the cor- 
responding chemical equilibrium constants, 

In the caiculutions in which constant composition 
flow is assumec to be maintained, the gas miature ree 
mains constant as the gases flow from combustion 
chamber to the exhaust atmosphere, Tie composition, 
thersfore, is equal to that at the adiabatic flame 
temperature (Tg) in the chamber, 

The chamber pressure (Pg) is taken to be 300 psia 

for all calculations presented here anc the nozzle 
exhaust pressure (Pe) is assumed to be 14,7 psia, 
Steady flow conditions prevail throughout the nozzle, 
It is agsumed that shock disturbunceées may be neglected 
and thet the velocity profile is flat, It is further 
assumea that one dimensional flow équations are 

valid, . 

Sinee the expansion process is taken to be isentropic, 
entropy chunges from viscous and friction effects 

are assumed to be negligible, 

The enthalpy of the reaction prouucts Is imiee 


pendent of the pressure, This is essentially 
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Se Neda kaoh weloolty (e*) io Gefinee arbi- 

j £6 Fesate Snes Slo), clamber |resture ans nod.de 

Wit Ores, Ad Such, 35 24 5 KeNMUFE of ted afficlency 

“of Clim comeustion process eas Can by experioortuliy etur- 

fem euodiy poumersa Quantitias, 4 ciworstionl vhlue 

Fw aes tin tv Se Sc, as Ty for » given propellant 
ust ceomistions, The thecretieul value is 

‘ Pelmsubit SF Ste propellant 4nd ia uswelly 

at ten percent Greater Cuen the eoperinent-liy detereined 

| = Comparison of the ebarscteristic veloeit, obtained 

ya os in tiwory, 15 4n eic in Combustion clanber 
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Pie effective exla 2 velovity (c) is tht velocity 
Wikich wien Mulbizile. by the made rate of Plow will give 
tie thrust, The effective exameust velocity te tle nocale exit 
Velocity heh conditiofe 4r@ a5 SssusGi in Uiie treats, Bere 





Ba 


again is 2 parameter which cun be computed theoretic.lly 
and also can be d@eternminea exzperimentully from easily 
measured quantities, 

Closely related to the effective exhaust velocity 
is the specific impulse (Igp). This is defined as the 
pounds of thrust per pouna of propellant consumec per second, 
It is obtained by dividing the effective exhaust velocity 
by the acceleration of gravity. 

It has been found useful in practice to define a 
nozzle thrust coefficient (Cp) which is a function of the 
' chamber: pressure, nozzle throat aPea, &nd the thrust, Thus, 
the thrust coefficient can be determined experinentully and 
compared with the theoretical value, The thceoretic.cl value 
is obtained by diviuing the efrective exhaust velocity by 
the characteristic velocity. 

The altitude index of a propellant system is the 
vertical range which a lurge rocket of an arbitrarily chosen 
. Gotal impulse would attain for a particular set of given | 
assumptions, Once the total impulse has been chosen and the 
assumptions made, the altitude index is « funetion of the 
‘specific impulse ana the mean propellant density, Therefore, 
the altituae index gives the relative value of a propellant 
system if used in a definite rocket which further satisfies 
the given conditions of manufacture and operation, These 
conditions assunle the use of a bipropellant syste, but the 
application to the systems studied here is still of value, The 
values of altitude index computed are only slightly higher 


than they should be, 
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DEVELOPULNT OF THE EQUATIONS YU SOLVING The CokPOO I TION 
OF THS PRODUCTS CF RSACTILAN INCLUJISG A SAMPL, CALIULATIO“ 


The method of wolving tlw composition of the pro- 
Sucts of réaction presented here 165 quite similar Vo thet 
uSed At tne Jet Propulsion Laboratory, However, it was er- 
rived at indejeraently sna it was useu for o1l coeposition 
solutions necessary in the préparetion of the results con- 
tained in tnie ChLUGLS. 





The propellant systees investiye ten contatives only 





three chemical vlementa; hydrogen, nitrogen, ata oxygen, 
‘The following SsyBbdols are exployed to represent the atcgic 
ana @olecular Specie. present in the Comoustion gases: 

= number of moles of water vapor (0) 

= number of moles of hyurogen (h,) 

= nusiber of moles ot hyuroxyl ion (OH) 

nusber of moles wf atomic hyuregern (hf) 


= number of moles of oxygen (0) 


~~ cC¢ A Oo CF BB 
a 


= nuwiber of moles of atomic oxygen (0) 
= nuweber of noiés of nitrous oxide (NO) 
® uumber of moles of nitrogen (Ny) 
nuzw ber of grap atoms of hytrogen 


* nusber ol gram etoes of okygen 


= oOo FF FF Ba 
8 


= nuwbéer of gram atous of nitrogen 
Thus, the gé@neral chemical equation for any systen 


Anvoiving only hycrogen, nitrogen, anu okygen becomes: 
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a 


m~ 


He WH + O = amyO + Biko ¢ COH + oh + CCe 4 TO + gi 
Lt can bu sevn that the probienm is to sclve for tle wight 
Ulbnown guantlebese a3 @, @, uy @y Ty By, ame By 

Prom this generai chemical equation, thret equations 


Lnvolvins «tom Dalances can be written: 


Sum of H atoms: He Za@¢ #b) €@ + (1) 
Sum of O atoms; O=aec+ze4f £ Bg Gap 
Suu of N atoms: R= gs 2h (3) 


The rewaining five equations necessary for — 
tion of the problem are obt.cined by bs tuming that complete 
equilibrius is reachea ana employing this faet to obtain 
five more independent relutionships among the eight unknowns, 
The equilibrium reactions involved anu the corresponding 


equations in terms of the unknown quantities are as follows: 


No + HoO = NO « Hp (Ang = 0,5) 
Kn; = gb | (4) 
hes 
«# H.0 = Op « cHe (Ang = 1) 
Kng * & p* (5) 
a” 
Heo = 0 4 He — (Any = 1) 
Kny © fb | (6) 
“a 
gHz = H (Ang = 0,5) 
kng = d (7) 
be? 
HoQ = OH + go (Anjo # 0,9) 


2” 
‘njo™ £2." (6) 
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In these equitions the equilibrium cometents are 
expresseu in terms of the mumber of moles of the compoients, 
In Tavle I the numeric. values of the evuilibrium cofistants 
are expressed in teres of partial pressures, K's, These 
were converteu to hp's by using the following general re- 


¢ 


lationship, 


Kn = K zl (6) 
where Dp = number of moles of prosucts of reaction 
P = total pressure of prouucts of reaetion (atz.) 
Qn = mumber of poles of ;roducts less number of soles 
of reactants as determines: from the balanced 
chemical eyustion for the particular eyquilibrius 
reaction 
Since the Kkyn's depend on the totsl number of poles of 
prouucts, we have in effect introduced a ninth unknown, Np» 


where 
Dp tasdbeceedesosfaegeh (10) 


which must «lso be estimeted for each particular solution, 

Assupbing a value for "np", these equetions can now 
be used to resolve the composition of the chawboer or exhaust 
gases, Rewriting (1), 


Caexbe-crwda2waHk 
From (68), ¢ = Knjo 8 
b° 
From (7), d ® Kno be? 
Substituting for "c"™ ana "d" in (1), 


2a + 2b s “19 a 4 Eng pro =H 
b* 


SBN a 


q™ 


i? 


lving above eGfuation tor "a", 


© 


a # H be? = kpg b - gbl.© (11) 











From (5), ¢ *= Enig (1) 
bed 
Fron (7), d= Kng be? | (15) 
From (6), @ = Kng a© (14) 
be 
From (6), f= Kn a (15) 
ameqeste: 
| b “a 
Frow (xz), g@O4c04te eT (16) 
From (3), h=z (N - g) (17) 
Repeating (10), Mp z= atrbsece+ ates fegeh (18) 
Repeating (4), Ens = pb (19) 
hes 


These nine equations, (11) through (19), and the 
general relationship, ‘Kquetion (5), were used to solve the 
composition at one'’hunureu degrev intervai: of tuemper.ture 
by using the following triel anu error procedure, 

1. Hstinate py" and cvolucte tic equiiLibriwy 
constants (hyjts) for the corresponsing chamber 
or exhaust pressure. Values of the equilibring 
constants are listed in Tubie I, 

2. bstimnate "b", anu using equations (11) through 
(17), soive for the other components, 

5S. From equation (16), fina "np", Compare this 


With estimated value of Step l, 


16+ 


4. From equation (10), Site hny:s COlpare this #ith 
the Value of Otep 1, The veluw of Byy of Step 1 
is in error only if "n." cotisete is incorrect, 

5. Repeat the process aging » new extinate of "nm," 
anc “b" until the caputec "n," agrees wth tie 
eatineted value, and bp, from equation (1¥) 
agrees with true value of step 1, 

6, AS & Check of NUiericel work, Eywetions (1), 
(2), amc (2) stoulu be uses, 

To Lllustrate this mwthou « suaple cuiculetion will 





presentes here, The System choown for the caleulation 








AS the tripropéllent systém consisting of 4 stoichiometric 
Bixcture of nitrogen tetrosiae «te hycrozine with one half 
molw of liquis hyGrojen ucuecd, The @gustion for this ree 
action at 00° K is: 
MOq * 2Mzlg + S.5Hy(1) © 4Hy0(1) & Og @ O.bHy (20) 
At an elévateu temperature; 
Og * 22H, @ O,ciz(1) © aleO(g) + Dey + COE 
+ GH + eOp + [0 + gWNO 4 hi_ 
From the left hand sice of (20), it ewn be sewn thet 
RH = 8,000 N= 6,000 0 = 4,000 | 
Tre prodiem is to wolve the componition of tie pre- 
duets of this récection When it occurs in a roteet Boter with ~} 
chamber pressure cf P, = ov psis at an estin tec Ciweder 
temperature of Te = Seuv® K. 
The results of tme solution ere prusefile. in Table A, 


In the first attempt .t sorution, "np" Was A¥stbe to be 


TABLE A 


ReouLty OF SakPLi. C~aLCULATION 





OF THu COMPOSLTION G? TH. PRODUCTS CF REACTION 


ror T = on 009 B; pe Pr. = Mel. ati, 





SOLUTION 2 2 
Step 1: Ny (est.) 7.960 7,960 
Kng 0.00607 0.00607 
Eng U,0027 Q,0Oz7d 
Kny 0.0074k 0.00742 
Kng O.17b5 0,175 
Knyo 0.0676 | 0.0678 
step «:! bd (est,) 0. Guu 6.750 
a a 5 4489 bss 
Cc O.<65 0,277 0.270 
a 0,157 0,152 0.163 
e €./0ee 0,062 0,057 
f @ 052 0,085 0,038 
g 0,110 OwOL1 0,064 
h 2.84 29Gb £4 066 
Ap 7,850 7.826 74068 
ptep 4: Kn. 0.0147 0, 003580 0,008LE 








= 
















7.88. am. "b™, 0,800, This solution gave « hng ( Sep 4) 
@uch greater then Kn, (step 1), indicating that tre -lret 
@stipate of "D" was too high, The m, (Step A) was lees 
“phan Np (Step 1). 

* The secom, attempt to votuin « correct solution Ws 
m ¢ With np *#.U.750. A New ebtinate of "np" wee not Bade 
to tne réletively lerge err of tie Piret solution plus 








the fact that tie suiutin Ae WG) Mere seeO bie to chinges 
ip “eV. ‘fkbs eolution gave « Lp, (step 4) less tle bn 
(Step 1), itdicating that "o¥ wis tuc low, again nm, (Step 8) 
mS less then np (Step 1). 

It wis seen fros the Cirst two .lvetpte at solution 





t ">" puct be between 0.600 am 0,70 anc "Hy* Buyt be 
between 7,650 «na 7,606. hence, "no" wes GssuneG to ve 
7,640 anu "b", 0,770, These vaives yieldeuw tue following 
results: 

Kng (otep 1) = »,00u0x Np (Step 1) = 7,640 
Kn, (Step 4) = O,0081z np (Step %) = 7.604 





AS & Check on the mechanics of the olution, suming 





fions of tis grup ateas ul M, 0, ae ff were Goce aici. the 
following result; 
Using @qu«tion (1), 
H@@aeebe @*e © = Biwi (om errer of .yul) 
Using eyuktion (a), | 
Oseec# xce4af se g = 6.00 (no error) 
USing equation (3), 


N= g-¢ wb = €,000 (no error) 







Hemcec, tre thir. solution, bakec Gh ah epetiood 


"no" of 7,620 uu "bY of U,.770, Was Shown to be the correct 


vOolution, uch COfmponwaut wus determine. tv witiin = G,002 


pty 


moles, No greater accuracy woula be warrunteu since errors 
a] 
innercnt in tie Valuedn Useu LOL tire beats of Lormatson anu 


£ 


for equilibriuse conmstunts are or thio oruer of magnitude 


lurger, 





-. (+= 


PAsT III 
SAMPLE CALCULATION OF PERFORMANCE PAKAMHLTERS 


In order to illustrate the methods used in celculat- 
ing the results presenteu in this thesis, the perfcrmunce 
parameters for the system consisting of a stoichiometric 
mixture of nitrogen tetroxide and erivenine with one-half 
mole of hydrogen suced are presented in detail, The equa- 
tion for this reaction at an elevated temperature is; 

No04 * 2 NoHge 0.5 Hy (1) ® 

a HeO (g) + Db Be 4 ¢c OR 9d H4 @ Op 
+f£O0O¢g NO + bh Ng 

The calculations are besea on tlw assumption that 
two gram molecular weights of hydrazine sre present, This 


gives u totel mass of rexctants, m, of 157,06 grams, 


Step A. Calculation of Chamber Tempersture (Tc) 


The first step in the calculation was to vetermine 
the adiabatic flame temperature in the echapber (T,), This 
equilibrium reaction temperature is that temperature for 
which the heat releaseu by the reaction (qi$) equals the 
enthalpy rise of the components of the reaction (} ny db Risto) 
Since both quantities depenu on the composition, the method 
used is as follows: 


a). Chamber temperature was estimated to be 3200° K, 


Chamber pressure is, is assumed, <0,41 ato, 


bd). The composition was solved for 3<00° K using 
the method outlinea in Part If, A Esp? and co were CoL- 
puted, These results appear in Table B, It is seen that 

a H5<00 is less than cs which showed that estimate of 

Te = 3z00° K is too low. 

cc). Second estimate was 500° K anu (b) was re- 
peated, The results for T, = 8300° K showed AHso0" 1s 
greater than ty ft is therefore between 400° K and 
5500° K. 

d). A linear interpretation of both Qoay and AH 
was used to find AT where 

Tg = 5200 « AT 
This double interpolation can be expressed as? 
AvNeet” ‘ AL a HE509 . 0 95200) co 


* O2Gay 5200 


- Qsy ) 
Substituting values from Table B, 
216,84 +4 Ab (<n8 «Sie ~ £16,846) = 224,75 
, Fy 
¢ AT(<10.70 = 224,75 
160 


AT = 34° K 

e). For use in subsequent calculations the composi-~ 
tion at the chamber temperature, > and A ison. were 
computed at this point. The composition at 5c54° K was ob- 


tained by linear interpolation of the compositions at $u00° K 
and at 33009 kK, 


atin 


TASLS B 


RESULTS OF UVAMPLE UsLCOLeTIon oF 
CRARBAN TEMPLAC TURE 


Tt (°K) &Z00 53500 3u84 
a 4,518 5,066 &.466 
b 0.770 0,685 0,726 
e 0.270 0,544 0,2¥6 
d 0,153 0.211 0.173 
e 0,057 0.076 0.064 
f 0,033 0.062 0.089 
z 0,064 0,004 0.071 
h &, 968 2,958 <,¥65 
Np 7,854 7,948 7,574 
A Hao 216.64 k cal 2e5,9% k cul als.o5 k cal 
gy 24,75 k cal ¥10.70 k cal 212,08 k eal 


where A 200 a y * (moles of compone.t) 4 = ( A E500 of component), 
‘es Q& (21) 
ay = Qr (products) - wy (reactants) (xe) 


T 
Table II list. values of 4 Haoq for each couponent 





versus temperature, 

Table III lists velues of neats of formation of 
the various components useu in evaluating ee of the chamber 
and exhaust gasés ana of the reactants, 





(c*) 





At this point the characteristic velocity, c*, was 
computed, This parameter is a measure of the combustion 
efficiency and is a function only of the propellant proper- 
ties, chamber combustion conditions, ana nozzle conditions 
from the chamber to the throat, 

By definition, 

CH « Po ft (23) 

: a 

By substitution in (<3), 

c# = ao 2 1 Yc & Ry Te (24) 

T(5q) Frag) Me 
By substituting the proper values of "g", "R.%, and 


the necessary conversion factors, (<4) beconess 


ct = Hoa 5 Se To ft seen (25) 
‘Sie Teo LO Mo 


if Te is expressed in degrees Kelvin anu M,, in grams/mol, 

In arriving at @quation (<4) from (235) only one as- 
sumption is necessary, The ratio of specific heat capacities 
of the products of reaction is assumed to be constant from 
the chamber to the nozzle throat and equul to 8@3 We being 
the ratio of the specific heat capacities of the products of 
reaction at the eyuilibrium chamber temperature, This is 
a reasonably good assumption, The pressure changes markedly 
from the chamber to the throat; but the temperature, anda 


consequently, & , Changes only slightly, 


az 6e 


In order to use (ed), Me, Ue, amc [~"(5,) were com- 
putea, The averege molecular weight of the prouucts of re- 
action at equilibrium chamber tei perature wes obt. ined 
first. The weignt of rwactents (m) wis Chosen to be 157,06 
grams enc the number of moles of chawber prouutte was com- 
putea in step A, 


i. = ; (23 
+. ) 
= . = Uo BO ad 

Me ' . i ¥40 go 1 


The isobaric heset cuapecity of the procucts of ree 
action wis found by consiuering tow enthilpy change of the 
products of reaction with the system in eyuilibdriue over 
the 100° K temperuture range nésrest Tg, This enthelpy 
change was then divided by tia product of the «versace nun- 
ber of moles present in tie cheaber and tue te@perature 
range (100° h) to give (Cp) os In this calcul.tion, 

(Gp), $3209 
np(100° x) 
where = § B3299 © (A KS9g = AlSo0) * (Gy = Say 

5<00 5500 


peo ea 


Substituting mumerical values from Table 1, Step A, 


(.7) 








¥ = (625 9% - 216,34) kh caél @ (2#4,.75 - ~10,.70) 
k cal 


@ 23.15 E cal 
2 (7,804 t 7,946) moles 
"P z 


= 7.891 poles 


wos] ow 


Using (#4), 


(Cy), » Reale x» i cai 
7,891 moles x 100° K 


29,512 cal moi~+ Onl 











mm 
¥o? (Cp), = (Cp), (<8) 
(Cy). (Cy), a Ry 
= 29,512 cal p72 Ox-t 
(29.512 - 1,986) cal mol 
Bo # 1.073 


The function (~! can be represented by the follows 
ing linear equation: 
t= ,1047 + ,0045 © (28) 
! = 54 t 
(= (Se) 01647 4 ,5048 (1,073) 


% ,6405 


Repeating (<5), 


e# = 246.01 4] So Te ft sect 
[~ '(&e) 1G Eig 
To 


2 4 5 TY ot 
= 946 01 1,0 x Onwo4 rt sec 
e 4006 10 x 15.948 








The calculation of the chamber temperature assumed 


that the reactants in the chamber compietely reactea and 
that they rétaineu in the chawber for a sufficient tine 


to permit the establishment of equilibrium concentrations 


at h— 


of all the components, When the guses levve tie ci ober 
ana flow through the noszle, there are two #uchenisse as- 
supeu possible, Hither the composition of the ges mixture 
remains constant as it flows through ti nozzle, or tne 
cowposition is changing continuously in accorcance with the 
pressure variations, temperature varistions, eno the core 
responuing eyjuillibrium comstants, In reelity, tome inter-~ 
mediate state is probably saintsained betweun tiese two extremes, 
In calculating Te it will first be assuweeud that equi- 
librium flow conditions exist in the nozzle «es well as in 
the chamber, At the end of the calculation it will be 
pointed out wnut modifications were applied to this method 
in order to compute the parameters assuming consteit composi- 
tion flow, 
Since isentropic flow is assumed, the exhaust 
temperature (Tg) Can be calculated from the following 


% « ve(pal ty ; ‘alt, (50) 


Te 
where tp 7 $ BT eal mol 7? O-2 (31) 
Mp (Te - Te) 

The correct Te can be culculated by the following — . 

trial ana error method, 
Te - 
&, Using this estimated Tg evaluate jhe, and Np. 
5. Use Equation (81) to evaluate Cp and then 


use Hquation (50) to calculete Tg. 


a2 


4, The correct solution is obtainea when the cal- 
culated value is equal to, or closely upproximates the esti- 
mated value, In the caloulations presented herein, the eal- 
culateu Te was accepted as correct when it was within ten 
degrees of the estimated value, 

The detailed computation of Te appears in Table C, 
The final result, Tg = Ke70 a was arrivea at by the fol- 
lowing procedure, | 

a) Exhaust temperature wus estimated to be 200° K, 
Using method outlined in Part Ii, the composition was soived 
for #c00° K end a total pressure of Pe = 1 atmosphere, 

b) With the composition known, At" ane ave. 


were cuiculated, 


_ . #200 SSF 


O04 
QOHzsp9 «42wwas computed in Step A. 


¢) The total enthalpy change of the products of 
reaction that eccurs between chamber and exhaust was conm-— 


puted using the following forrula: 


3234 S234 2200 5254 
& Ho300 = ( HES 5 ~ AH359 ) + A®avozog 


d) The average number of moles of reaction products 
that are present aus the gas flows from the chanber to the. 
exit was found, 

ip - nowo4 t new90 

a 
e) Equations (30) and (51) were then used to com- 


pute Cp and Te. 


“ 


TABL® C 


AOBULTS OF SAMPLE CALCUL. TION OF 
EAHAUST TRMPERATUNE 


T(°K) 3434 (To) E=00 =300 
a 3.468 & WBS 3,969 
b 0. 7¥e 0,503 0.510 
C 0,295 0,012 0,04 
a 0.173 0.011 0,019 
@ 0,064 0,000 0.00% 
f 0,039 0,000 0,002 
g 0,072 0,003 0,002 
h 2,865 2,998 “999 

Op 7,874 7,512 7,686 
b B99 219,83 1Se,%b 140.40 
d Wve 9,00 44,84 43,44 
i 2 3 
2200 2u80 2270 
~19.85 219,85 218,85 
132,265 137.14 157,96 
44,94 44,04 43,89 
132.62 16.83 1U5,86 
7.693 7,697 7.697 
1034 974 964 
16,675 16.917 16 , 862 
2257 2270 £270 








Repecting (61), 


= 5234 
Cp = S$ H299 cal mol) x~1 





“2 } 
7,698 moles x (S254 - 2200) OK 
-1 O.-4 


Cp = 16.673 cal nol” 
And using (30), 


R 
fe = Te Pe| tp 
Fe 


1.986 
2 254° Kk 1 atm i 0. 


20.41 atm 
Te = eeb7° K 
f) The next estimate of exhaust tewper:ture was 
2260° K, Since this estimate was not am even hundred ae- 
gree value, the following were calculated prior to repeat- 
gy (b), (ec), (4d), and (e): Composition at 2500° K, 


ites Gn hE 
Akzog » and XQave300° 


g) The procedure outlinea in (b), (c), (da) and (e) 
was repeated for the new estimate of Te = 2260° K, Linear 
interpolation between 2200° K ana 500°K values was used in 


determining b HgRBO ; AQayocce) and . 


h) This second attempt at solution yielded Te = wet 
The variation from the estimated value was on the border- 
line of the accuracy standards set up for the celculetions 
so a third solution was wade using an estimated Te @ 2270" K, 
This gave a calculated exhaust tenperature of the same velue, 


Hence, the correct Te was taken to be «270° K, 


It shoula be notea at this point that 5 tig’ as 
required for the calculation of the effective exhaust ve- 
locity and pust be culculated using the exact uxlhmwust taepera- 
ture, not using u pengera sure witrin tem degrees of tre cor- 
rect exhaust temperature, fherefore, wien the caicul«ted 
and estimatea exhaust temperatures were not exactly the same 

 Walue, a reevaluation of 53° was made using the correct 
value of the exhaust temperature for use in evalusting the 


. effective exhaust velocity, 


(e) 








(Igy) 


The effective exhaust velocity was computed from 
the following formula; 
(32) 





where J is the mechanie:.] eqguivuiient of heat, 
In obtaining Bquation (32) it was ussuser thats 
1) The products of reaction behave as perfect gas, 
&) The expansion precess in the nozzle is adiabstic, 
5) The velocity of the gases in tie chamber 15 so 
small in comparison with "c" thst it can be 
set equal to zero, 
By substituting the proper vulue of J, anc the neces- 
sary conversion feetors (S5&) becomes: 


5 HAC ft see) (33) 
b 







@ ws 9495. 


~632 


Te 
if 6 Hy is expressed in kilocelories anda "m" in grams, 


rT 
The § Hip, was computed in step C and mw wus chosen 
arbitrarily at the beginning oi the celeulation, Substi- 
tution in (63) gives; 


¢ = 9495.2 1:35, 86 Lt see~+ 


157.07 
¢ = 8457 £t sec7t 
The specific impuise is expressed by the following 
relationship: 
Isp = ¢ : (54) 


w 8497 ft sec? 
33.2 £tb see” 





Igp * ®605,9 sec 





A nozzle thrust coefficient is aefined by the fol- 
lowing relationship: 


rs _£... (35) 


Expressing tue effective exnaust velocity in terms 
of the thrust and mass flow rate, it is seen thet 
Caw Ff =z F/Pe¢ fy = Cr 
~. ak 3 eit 
m u/Pe ft Cz 
or 


Cr = ¢. (56) 
c+ 





~O4- 


ft seem) 
ft see~1 


oy = aay 


Cp = 1.39 





The average moleculur weight of the reaction pro- 
ducts curing their passege through the nozzle was computed 
from the relationship, 


T= z (37) 
where Tp = 3254 ¢ 2270 


2-2 
- a tie. moles 


Wp = 7.627 moles 


Ws golad on — 


= £0,407 gb pol~t 








The determination of altitude index was made from 
a set of curves which weligntec the specific impulse ani meen 


propellant density in a mainer believed more properly than 





given by déensity-impulse, This altitude index is for & 
large rocket manufactured and operated under a given set of 
assumea conuitions, 


The mean propellant tensity (+ ) was calculatea, 
' VoTume 
; 
poe a 


- = 1,lo gm on@9 





aS = 


The values of the densities used to determine the 
totul volume occuplea by tuls mass of Liguia propellants 
are listed in Table IV. 

From 6tep D, Igp s «00,0 Sec, Using this Igp and 
the ~ calculutea ebove «es entering argusents to the curve, , 
the ultitude index was founu to be 556 miles, 

For the case of constant composition flow, the 
chamber temperature (Tg) and the characteristic velocity (c*) 
are Computed exactly 4s outlined in Steps A anda B, since 
these quantities are independent of the type of flow in 
the nozzle, 

In calculating the exhaust temperature (Te) in 
Step C, the composition of the gis mixture was held fixed 
at the composition that pPevailed at the equilibrium chamber 
temperature ana the changes in the resulting heat balances 
made, As a result of this constant composition, d avn 
Will vanish and § Hg’ ’ Cp and Te will be lower than the 
Values obtained assuming equilibrium flow, 

With these exceptions the remaining calculations 
are similer to those for the case of euyuilibriua flow, 
However, the values obtaineu in Steps D, E, ana F will 
aiffer from those obtained in the equilibrium flow caleu- 
lations since they uepend on 5 Hee and Np Calculated in 
Step C, The characteristic exhaust velocity (c), specific 
impulse (Igp), nozzle thrust coefficient (Cp), «nd the al- 
titude index (h) are ail less than the corresponaing values 


obtuinedad assuming equilibrium flow, The average molecular 





weight of the rusetion products (I) is grester than the 





St 
: PART LY 
DISCUSSION OF RESULTS 


The performance parameters cf the three bipropellant 
systens in stoichiowetric proportions, which were usea as 
the basis of the tripropellent systems stuaie. in this 
thesis, were known earlier (Cf. Refs, ©, 11, anu 12) al- 
though corrections were made to the hydrogen peroxide- 
hydrazine system, These corrections were made since 
earlier calculations of the performance parameters of the 
hydrogen peroxide-hydrazine system (Cf, Ref, 9) neglected 
molecular and atomic oxygen as components of tie combustion 
gas mixture, All three systems investigatea exhibit rela- 
tively high performance and, as such, they held promise | 
for future development and are currently of interest, 

The problem of aetermining whet effect the addition 
of drogen or ammonia woulu have on each of these systems 
wus investigated by culculating the perrornuance paraneters 
for two series of tripropellant systems, Tie first series 
consisted of the stoichionetric mixture of the bipropeilant 
systen, to which was added increasing amounts of hydrogen}; 
the second series, the stoichiometric mixture of the bipro- 
pellant system, to which was added increasing amounts of ume 
monia, for each systen all the pursmeters illustr.ted in 
Part III were culcul:ted for both equilibrium ana constant 
composition flow conditions, The results are tabulczted in 


Tables V~XVI and are presented graphically in Figures 1-16, 
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For eeck tripropeliant system studleu tie Jollow- 
ing results &re presentéu: 4 table of all psraheters Com- 
puted assuming equilibrium flow comitions, a table of 
all pseraweters computed assuming constant composition flow 
conditions, a ¢ruph of chasber te@perature, exhaust tempera- 
ture, specific igpulse, anu Gititudé index versus weight 
per cent of thiru component for both equillorium ami con- 
Stant composition flows, am a graph of the percentage 
variations of these sam@ paraoeters with weig)lt per cent 
of thira component, 

The graphs of the percentige veriations of certain 
of the purupeters were prepsrea in order to better illus- 
trate the relative effect of the thira camponent on the 
Various parameters, These graplh# are based on the ratio 
of the parameter for s given percertage aduitive to the 
value of the parameter for the stoichiovetric mixture of 
the bipropellant system, by referring to these curves for 
any particular system one can tell what the effect on the 
chamber temperature, the exhaust temperature, specific in- 
pulse ana altitude index will be for a given weight per 
cent of third component - all in terms of a perceiitage of 
the value of the paraseter before the addition of a third 
component, 

At the first inspection of « plot of percentage 
Variation of paraseters with weight per cent of third com- 
ponent, it may appear that the constant couposition velue 


of the parameter is greater than the corresponcing equilibrium 


value, It should be realized thit sithoug:. the percentage 
of the stoichiometric value uay be greater for constant 
composition flow, the absolute velue of thé paraneter for 
equilibrium flow is siways greater tian, or at least equal 
to, the value for constant composition flow, 

A, ADDITION OF LIZUID H¥YDROG:LN TO A STOICHIOKLTRIC 


KiaLUk. OG. NIVROGEN TLTROALDiC AND HY¥DRacIni., 


Of the three basic bipropelilant systems the best 
performance is shown by the nitrogen tetroxide-jhydrazine 
systen (Cf, Ref. 11). Its chanber temperature at the stoi- 
ehionmetric mixture ratio 1s 5255° kK, Assuming equilibrium 
flow conditions the exhaust temperature is 2334° K, the 
specific impulse is <o9.2 seconds, and the altitude index 
is 5<5 miles, Corresponding values for constant composi- 
tion flow are 18<5° kK, #47,5 seconds, and 47< miles, respec= 
tively, These values are for the parameters at the stoi- 
chiometric mixture ratio and are tne basis for the compari- 
sons presented in succeeding paragraphs, 

The chamber yn exhaust temperature,: spe~ 
cific impulse, and ultitude index for both equilibrium and 
constant composition flow conaitions versus the weight per 
cent of liquid hydrogen in tne tripropellant system are pre- 
sented in Figure 1, Figure ¢ is a graphical presentation 
of the percentage variation of the same parameters as 


plottea in Figure 1 versus weight per cent of liquid hydrogen, 


= 


=4\ pe 


The values of @ll parameters computed as well as the per- 
centage Variations plotted in Figure « are listed in Table V 
for equilibrium flow and in Table VI for constant composition 
flow, 

Figure 2 shows that the greatest effects of the hy- 
drogen addition are sewn in the changes of the exhaust 
temperature anu the Chamber temperature, Por ejuillbrium 
flow the exhaust temperseture decreases steadily with ine 
ereasing hydrogen, For tue systee containing paxieuw hy- 
drogen, 16.25%, an exnaust temperature of 509° K is only 
54.7% of the stoicihiowetric value, U3549 }, Tre effect 
On tné exhaust temperature for constant composition :low is 
hot quite as great. At the point of Maxigus hydrogen adci- 
tion the exhaust tempersture of 809° & is 44.4% of the 
stoichiometric value, lseé° kh. The echseber temperature 
varies with hydrogen addition in the same cenersl] Banner as 
the exhaust temperature, The effect, however, is not as 
pronouncéd as it is on the eXmaust temperature. The @axizus 
effect is at the point of the greatest hydrogen adiition 
where the chamber temperature is still 1692° K, 52.3% of 
the stoichiom@etric value, S3258° K, 

Thus, it can be seen thet the effect on the chagber 
ano exhaust temperatures 18 very desireable ano appreciable, 
The next consicerations are the corresponiiing effects on 
the performance parameters, [t can be sewn from Figure «& 
that the effect on the specific ippulse is also cesirable, 


Tre addition of hydrogen increases tme specific ippulse 
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which reaches its maximum near the point corresponding to 
9.0% hydrogen, For equilibrium flow the maximum is 275,86 
seconds, 106.5% of the equilibrium stoichiometric specifie 
impulse, 269,2 seconds, For constant composition flow 

the maxinum is <75,.4 seconds, 111.0% of the constant com- 
position stoichiometric value, <47,.5 seconds, Both curves 
are very flat, varying less than 4,5% in the range of hy- 
drogen from 5 to 16%, 

The effect on the altituce index is comewhat dif- 
ferent, again referring to Figure <, it is seen that for 
low values of hycrogen the initial effect is to increase 
the altitude incex, The maximum is soon reached, and then 
the altitude index decreases in much the Same manner as 
the temperatures because of a rapidly declining mean bulk 
density, ‘For equilibrium flow the maximum, 556 miles, 
102.1% of the stoichionuetric value, 5<5 miles, is. reached 
at 0,64% hydrogen, The constant composition flow maximum, 
496 miles, 105.1% of the stoichiometric altitude index, 
4%< miles, is reached at <,5% hydrogen, From these maxi- 
mum values the altitude index decreases steadily anu for 
the system containing 16,<35% hydrogen, the values are 
580 niles, 62.9% of the stoichiometric value for equilibrium 
flow and $36 miles, 70.0% of the stoichiometric value for 
constant composition fiow, 

in spite of the decrease in altitude index, the 
overall effect is still desirable since the percentage 


decrease in chamber temperature is greater than the corres- 
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ponding percentege Gecrease in sititude index, In prectice 
the two faetors woula have to be carefuliy weightec depernm- 
ing on the use to which the missile is to be put. For hy- 
dregen aduition greater than 3), the percentage decreose in 
altitude index lags the percentage decrease in chatber 
temperature by 4 to 10% for equilibrium flow; for constant 
composition flow, by 1l to 1e%, 

The characteristic velocity varies in ualrost the 
exact manner in which the effective exhaust velocity and 
specific impulse vary. This can be se@n by examining the 
variation in thrust coefficient, the ratio of effective ex- 
haust velocity to the characteristic velocity. Por equi- 
librium flow it varies from 1.40 at stoichiometric to 1,36 
for maxigums hydrogen. This is 1tss than <% varisetion, For 
constant composition the veri«tion is sligi tly greeter, 


from 1,54 to 1.58, or about S%, 


B, ADDITION Of LIQUIL ARWONIA TO A STOICHIOME,TRIC 
MIATURE OF NITROGEN TETHOXIDez AND HYDRAZINE, 


The results for the syste, nitrogen tetroxice- 
hydrazine with liquid ammonia auced, are presented in Table 
VII, Table VIII, Figure 3, anu Figure 4, For ;urposes of 
discussing this system, Figure 4, the percentage variation 
curve, will be used, 

As in the previous system discussed, the greatest 


effects of the third component are on the temperstures, 


we ae 


The chamber temperature ana exhaust temperatures uecrecse 
steadily with increasing ammonia, At the point of maximum 
ammonia addition, corresponding to 17.91% by weight, the 
temperatures are: chamber temperature, <626° K, 61.5% of 
the stoichiometric value, 3232° Kk; exhaust temperature for 
equilibrium flow, 1456° k, 62,4% of tie stoichiometric 
value, 25549 Kj and exhaust temperature for constant con- 
position flow, 1434° K, 78.7% of the stoichioxetric value, 
1625° K, 

Tne specifie impulse, assuming equilibrium flow, 
decreases slowly with increasing ammonia ana for the maxi-~ 
mum ammonia content, the specific impulse is 244.6 seconds, 
94.4% of the stoichiometric value, 259,2 seconds, For cone~ 
stant composition flow the specific impulse increases to a 
maximus of 250,8 seconds, 101,5% of the stoichionetric 
value, <47,5 seconds, near 5% ammonia, It then decreases 
to <42.7 seconds, 98.1% of the stoichiometric value,at 17,914 
ammonia, 

Similar sendecaiite are shown in the variation of 
the altitude index, For equilibrium flow the altitude index 
is almost constant from 0 to 4% ammonia anu then dreps off 
to 439 miles, 85.6% of the stoichiometric value, 525 miles, 
at maximum ammonia content, For constant composition flow 
the altitude index increases to a maximum of 4&5 miles, 
102.5% of the stoichiometric value, 47 wiles, near 5, 
ammonia, It then decrezses to 44 miles, be,6% of the 


stoichiometric, at the point of maximum ammonia, 
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As in the system in which hydrogen was added, the 
variation in the thrust coefficient is very suxll indicating 
that the characteristic velocity varies in almost the same 


manner as the specific impulse, 


C, GENERAL RESULTS APPLICABLE TO ALL SYSTERS 
INVESTIGATED, 


Two systems have been discussed in detail, Similar 
discussions of the other four systems studied would not be 
very dissimilar although numerical results are not identical, 
The results of all six systems are presented in Tables V 
to XVI and in Figures 1 to 16. A study of these results 
makes poSsible certain generalizations to all the systenis 
investigated, 

When the per cent variation of chamber temperature, 
exhaust temperature, specific impulse ani altitude index 
are plotted versus weight per cent of third component, hy- 
drogen or emmonia, the following can be said of these curves 
for all systens considered, 

1. TLe parameter that is affected to the greatest 
extent is the exhaust temperature for eyuilibrium 
flow, 

2, With the exception of the red fuming nitric 
acid-anmonia-hydrazine systen, the parameter that 
is affected to the next greatest extent is the 
exhaust temperature for constant composition 


flow, 
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With the exce,tion of the red fuming nitric 
aCid=Mmonia-hyurazine system the paraneter 
that is aftected to the neat greatest extent 
is the chuuber temperature, Both 2 and 3 are 
also true for ammonia aduition of greater 

then 6% to the red fwming nitric acid-a.ronia- 
hydrazine systeus, 

Next in order is the altitude index for eqgui- 
librium flow, However, there are several 
minor irregularities in this order for weight 
our Saute of third component of less than 2%, 
If the altitude index for constant composition 
flow is neglected for tue time being, the 
parameter which maintains the highest per~ 
centage of its stoichiometrie value is the 
specific impulse for constant composition flow, 
Just below this is the specific impulse for 
equilibrium flow, 

fhe altitude index for censtant composition 
flow belaves quite differently with the addi« 
tion of liguia hydrogen or liquid ammonia, 
This occurs because of the large difference 

in the densities of the two liquids, Liquid 
hydrogen has a density of .07 grams per cubic 
centimeter while liquid ammonia has a density 
of .648 grans per cubic centimeter, Por 


hydrogen addition of greater than 5%, the 
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altituae imdee for constant composition f low 

is below tite curve of the specific lmpulse 
variation. For ammorls. aduition this paraneter 
remains above the Specific impulse curve for 
constant composition flow out to a point where 
the ammonia etdition is pdetween 6,5% and lig, 
depending on the system, 

The main purpose of the aWuition of & thiru compo- 
nént is to reauce the chehber Lemperature, Howéver, with 
this change of Cemperaturée une the presence of incrwssi'g quan= 
tities of « third component, the performance clanges, Ip an 


attempt to determine parawetery which can be use. 485 & Bedsure 





of the effectiveness of the wauition of the thirc component, 





the percentage variation curves were considered, eno frow 
them the following quantities were computes for “all systenas 
for both eyuilibrium and constant composition flows: 

E(Isp) defined as (Isp/Isp° - Te/Te°) 

E(n) defined as (b/h® - Te/Te-) 

The superscript zero represents the value of the 
parameter for the stoichiowetric mixture ratio before the 
addition of a thira component. It was found that the value 


of each of these paraweters was practically constant for 





equal percentages of hyurogen in the three systems, The 
same was found to be true for eyual percentages of ammonia 
for the ammonia containing systeis, 

These values were averageu for tie hyvrogsn systeus 


énu for the ammonia systems, The average error in the velue 


or i I er B Ch! for any syStem, «en competed with bie 
Sp i) 
- t 4 


mean Vaiwe of the Purkaedter Vor tle une 29Sttr:, war 


found te be small as tan de Been Trem the Tollowing resuies 


for the nitrogen tetroxide~hydregen-hydrazine system, 


Elow Paraneter AV. Error Max, brror 
Equilibrium “(Isp) 0,004 3.008 
Constant Composition E (Isp) QO.015 0.021 
Equilibrium Eh) 0,006 0.007 
Constant Composition E(h) 0.087 0,050. 


This order of magnitude of the «average anu maxitmm errors 
wa, the seme for the other systems, 

The velues of these perumeters are listed in 
Tables AVII to XXIV, Figure 17 is a plot of the averuge 
vaiues of the four quantities, (Isp) and “(h) for both 
types of flioew, versus weig.t per cent of hydrogen or am-= 
monia, A study of Tables XVII to XXIV and Figure 17 makes 
several general results apparent, 

For equal weight per cents cf hydrogen, the jer 
cent lag of the specific impulse or altituace index behind 
the per cent variation in chamber tempéruture is the sune 
within # 3% for all three hydrogen tripropellant systens 
studied. This is also true for equal pereentages of ammonia 
in the three anmonia tripropellant systexs investigated, 

The obvious usefulness of these relatively simple 
parameters lies, of course, in their ability to predict 


performances of tripropellant systems onee the chamber 
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temperature curves have been c.iculuted or estineted, 
and a5 an indirect means for iliustrating the relotive 
merit of « series of thiru components uncer consiuer.«tlon 
as possible coolants, 

The parameter, E(Igp)s WGS Qppliec to uste frow 
Other sources in so far aS it wes aveileble, I[t is shown 
in Table XVIII that &(Igp)? Calculateu frow data for constant 
composition flow conditions for the liquid oxygen-liquid 
hydrogen-hyarazine system (Cf. wef, 3), wes in good agree- 
pent with these results, The average error was ,016., This 
was the only other data available for « tripropellent system 
containing hydrogen, nitrogen ana oxzygen., 

Using dats for carboneceous triprojpellant systeses 
(cr. Ref, 4), it was found that agreement among the car- 
bonaceous systems was satisfactory, but that E(I,,) for 
systems containing carbon was not in agreement with E(Isp) 
for those tripropellant systems containing only hy) crogen, 
nitrogen anu oxygen, However, it is believed tiat sitilar 
parameters might be useful in analyting carbonaceous sys- 
terms at such tiwe as more tripropellant data is «veiluble, 

The relative effectiveness of hydrogerm enn be con- 
veniently comparéd with that of ameonis on th® basis of 
these parameters, If the varietion of specific impulse 
is of primary importance, “(n,,) is the interesting psraneter, 
For low pereentage of hyurogen (say 4%) E(I,p) is about 
nine times greater than the correspon ing “ (Isp) for ss@onila, 


At highér percentages the etrect is less pronouncec, For 
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example, at 15% of hydrogen the L(y oe is about four tines 
that for anmuonia, é 
When the variution of altitude index is of primary 
importance as in the case of long range or high altitude 
missiles, hydrogen is also more effective than anmuonia + 
based on a comparison of corresponulng Values of E(h)- 
For 4% of hydrogen or ai.monia the ratio of the E(h) for 
hydrogen to the Ein) for ammonia is eight, For 15% of 


third component the factor is still two, 


10. 
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A-3S 
TABLE IIL 
HiaTS OF FORMATION USED IN THis INViscTIGaTION 


Ho(L) eeccsecees 24648* k cal noi~t 


Hgbg wtacavisees S808" ® 
HoOg(4) voccnwe. 328.2" r" 
HNOg sessccoeeee #41,66° rn 
RFNA(6 8% 
NOg eevee #61,00 " 
NivHg ceoseccvesee “12,05 ft 
NH (1) sewesecee 416.07° fn 
HoO(B) wesseeees #576798 i 
GH sssisnssscsass @LOS08 nt 
DO cesesesecesese “09,109 i 
WM uct cbbesesstis Meke ORE i 
WO ievssewueswen thbeOF mn 


fcalculated by L. G. Cole from data in Chemical Rubber Pub- 
lishing Company Hunabook of Cheriistry and Physies, pp, 1745- 
1747 (1945). 


Dchemieal Rubber Publishing Company Handbook of Chemistry 
and Physics (1947). 


“Bichowski, F, R., Rossini, ¥,. bD., "The Thermochenistry of 
Chemical Substances" Reinhold Publishing Corporation (1936), 


All other data from "Tables of Selected Values of Chewical 
Thermodynamic Properties", National Bureau of Standzrds 

(U. S. Department of Commerce) with support of Office of 
hiaval Research, USN, March 81, 1947, 
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TAILLE IV 


DENGITIBS USED IR THIS INVUSTIGATION 


Component 


Wioig 

HO, 

NOg (1) 

Hy (2) 

NH» (1) 

HNOg (6.8% NyOg) 


Tr 

gm con” 
1.01 
1.465 
1,41 
0.07 
0,648 
1.548 


at 18° ¢ 

at o%¢ 

at o° ¢ 

at «BbE,6° C 

at 20° c (Cf. Ref, 14) 
at 20° c 


b-5 
TABLE V 
VanlaTION OF PRRFORM.NC: PaARAMATERS WITH WEIGHT PER CONT 
OF LIQUID HYDACGEN FOR NyOQg-Hy(1)~NeHa SYSTEM 


ASSURING EQUILIBRIUM PLOW CONDITIONS 


(2) 









(i) (S) (4) (5) (6) (7) 
HES HOLES rT. n/i” | Oo Te/Tae Ff 
°K 7K gm em 
0.00 0,00 3255 1,000 2334 1,000 1.25 21,446 
0.64 0,50 6e84 1,000 2270 0.975 1,13 20,407 
1.284 1,00 $208 0,992 4149 O.9%l 1,08 19,418 
1.80 1,50 S108 0.977 <O057 0.875 0,94 186,452 
2<s2 te Sage 0.969 i930 0,827 0.66 17,597 
4,91 4,00 2809 0,868 15u8 0,685 0,68 14,866 
8.37 8.00 “285 0,706 1185 0,605 0.48 11,475 
25 W100 i6S2 0,53 809 0,347 0,83 8.468 
(3) (10) (12) (le) = (18) (4) (28) (26) 
WT Cr e Isp Isp/Isp® Cp hb h/h?- 
ft sect it see7+ sec mi, 
0,380 oGGH 8546 <bu.2 1,000 1.40 625 1,000 
0,64 6108 _ 8497 {65.9 1,018 1.59 6556 1.021 
1.28 6198 8587 266.7 1,029 1,38 527 1,004 
1.90 6264 86bb 266.6 1,087 1.36 521 0,992 
se bie 6508 871s 270.6 1,044 1.48 518 0,989 
4,91 6590 68a 2741 =, 087 1,38 492 0,987 
9.57 6413 3u68 <7b,4 1,065 1,38 4x0 0,860 
16.25 6557 UTS S71.3 96-4 087 1.58 330 0,629 


VARLTAZION OF Pah WWANC 


(1) 
Ws 


0.00 
0.64 
1.<8 
1.0 
& Se 
4.91 
9,57 
16.235 


(9) 


0.00 
0,64 
1.28 
1.90 
240% 
4.91 
9.57 
16,235 


(2) 
WOLES 


0.00 
0.50 
1.00 
1.50 
«00 
“00 
8.00 
15,00 


(10) 


Cc# 


0965 
0100 
6198 
6x04 


Anti 


TABLe VI 


(4) 
t./Te 


1.000 
1.200 
0. Bex 
9.977 
0,859 
0.668 
U. 706 
O.,b25 


0 


ASSUMING CONST NT COMPOSITION FiO! 


(5) 
Te 
°K 

16e5 

1621 

1601 

1768 

1727 

15.5 

117% 
60S 


(6) 


,/t," 


1.900 
0,09 
0.956 
0.“70 
O,u47 


' 0,657 


),645 
0,444 


OF LiylUly HLUnOGEN FOR NpOg-Hp(1)-8.H,g OYSTEE 
COWL. TIONS 


(7) 
e 


1.9% 
1,16 
1.05 
0.04 
0.68 
v.68 
» 468 
0,35 


gr en7> 


PAAARUTLAS GITH OVLGNT PAR CENT 


(@) 
7 


<0, 200 
i848 
1#.117 
16.005 
17.410 
14.611 
11,470 
&.4665 


ft secm+ ft sec-t 


7870 
8161 
65<4 
64.56 
GSd2 
B76z 
6867 
6757 


(12) 
Isp 
sec 

247.5 

253.4 

258.5 

wUw .0 

269.0 

Teal 

270.4 

<71,.5 


(15) 


Isp/Isp° 


1,000 
1,04 
1,044 
1,009 
1.071 
1,099 
1.115 
1,036 


(14) 
Cr 


1,54 
1,34 
1,54 
1,55 
1.35 
1.57 
1.56 
1,08 


(15) (16) 


h 
mi, 
472 
468 
491] 
4906 
496 
481 
420 
350 


h/n® 


1.000 
1,036 
1.040 
1.061 
1,051 
1,019 
0,820 
0.700 
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TABLE VII 
VARIATION OF Pi.RFORMANCE PARAMLTERS WITH WEIGHT PER CENT 
OF LIQUID AMMONIA FCR NgOg-NHg(1)-N2H, SYSTEM 
ASSUMING EQILIBRIUK FLOW CONDITIONS 


(1) (2) (3) (4) (5) (6) (7) (8) 
"% Mb tf iD & /%"” 7 ¥ 
OK OK Zr em? 

0.00 0,00 3283 1.000 #334 1,000 1,25 23,446 

2,66 0,25 5208 0,998 2246 0.962 1,22 20,842 

5.17 0,50 161 0.978 2102 0.901 1.19 20,214 

9.83 1,00 S017 0,933 1842 0.789 lyl4 18,071 
14,04 1,50 #2830 0,875 1628 0.698 1.10 18,088 
17.91 2.00 <626 0,813 1456 0.624 1,07 197,255 

(3) (10) (Li) (12) (15) (14) (15) (16) 

wT ct C Isp - Isp/Isp° Cp h h/h® 

rt see"+ ft seec+ sec mi. 

0.00 5965 8546 2b9.2 1,000 1.40 55 1,000 

2,66 6014 e375 260.1 1,008 1,88 527 1.004 
5.17 ~~» «6017 8343 259.1 0,999 1.39 523 0,996 

9.83 5976 8229 255.6 0,986 1.38 498 0,949 
14.04 5848 6073 250.7 0.967 1.38 472 0,899 
17.91 8705 7875 244.6 0,944 1,38 439 0,836 


£8 
TABLL VIII 
VARIATION OF PLKVOMMANCE PARAMETERS WITH WHiGHT PLA CENT 
OF LIQUID AMMONLA POK NoOg-WHg(1)-NoHg SY¥UTeD 


AUSUMING CONSTANT COMPOSITION FLOW CONDITIONS 
(1) (2) (3) (4) (5) (6) (7) (8) 
W™%% $#MOLES = § Te fi.” & /t,.° Ff ¥ 
°K °K ge om~* 
Q.uV0 0.00 5255 1.00 18xv35 1.Q00 1,25 £0,800 
2.66 0.¥d Szu08 0.982 1606 0,90 1,be £0,564 
5.17 Y.50 5161 0.¥76 1760 0.976 1.16 19,886 
6.63 1.00 $017 O.955 1687 O.Beb 1.14 16,81) 
14,04 1,50 W550 0 «0 B7HCiéidSSC(‘i‘i ‘dw; OCOL10)—s-16,0918 
17.91 2.40 ~6<8 0,815 1454 0. 787 1,07 17 .<k3 
(9) (10) (11) (lz) (18) (14) (15) (16) 
WTs ce Cc Isp I, p/ I ~ Cr h h/h® 
ft secm+ ft sec-) sec wi, 
0,00 5965 7370 247.5 1.000 1.34 472 1.000 
<.66 6014 6059 249.6 1.008 1.4 47S 1,015 
5,17 6017 8076 250.8 1,013 1.54 483 1,028 
9.835 5976 8065 250.5 1.012 1.35 €79 1,015 
14,04 9648 7286 243.5 1.008 1.37 460 0,975 
17.91 S705 7616 24¢.7 0.981 1.37 44 0,698 


A-@ 
TaBLbs ITZ 
VARIATION OF PARFORMANCK PARAnseTARs (1TH WetGhT PBR CENT 
OF LIQUID HYDsOGEN FORK H,Ou-Ho(1)-NHg SYSTEb 
ASSUMING & ,UILIBRIUL FLOW CONDITIONS 


(1). ) (5) (4) (9) (6) (7) (8) 


Whe MOLLE To [i . > at OF Y 
Oy °%K gm em7° 
9.00 0,00 6bl 1.000 15% 1,000 1.28 ip 7a? 
1,00 0.50 <845 0.997 1605 6,912 1,08 186,250 
1.97 1.00 2748 0,964 1670 0.844 0.55 16,951 
2.938 1.50 +639 O,9%6 16844 0.780 0,85 15,8351 
3.87 2.00 2527 0,686 1445 0.730 0.77 14,654 
7.46 4,00 <140 0.751 1244 0.578 0,56 12,013 
13.88 5,00 1626 0.571 788 0.403 0,58 8,938 
25,21 15.00 1145 0.401 588 0,267 0.26 £6,515 
(2) (10) (11) (1) (15) (14) (15) (16) 
Ly cw c Isp Isp/Isp”? Cp bh  h/h® 
rt seewt ft secm sec mi, 

0,00 SITS 3067 280.5 1,000 1.40 481 1,000 
1.90 LOSE Bic 266.5 1,008 1.89 487 0,992 
1.97 95976 BzeO -«eb7.5- 1,028 «= «SS 47S (0, 867 
2,9 6012 B54y £69.3 1.045 1,39 462 0,941 
5.67 6044 8546 260.7 1.041 1.39 456 0,929 
7.46° 6095 6461 i8<,8 1.049 1,39 409 0,833 
15.88 6085 8596 260.7 1.042 1.88 S26 0.664 
25.u1 5886 8128 0be.4 1,008 1,88 227 


0,462 


A-10 
TASLL X 
VaR aTIO“N OF PLAPWRAANCE PARAMETERS WITH WEIGHT PER CERT 
OF LIQUID WYDROGLN /On HyOy-Hy(1)-ByR, SYSTEM 
aAVeUNING CONSTANT LURPOSITION FLOW COMUITIONY 


(1) (<) (S) (4) (5) (6) (7) (8) 
wT MOLES =f, ./%" & G/" F Y 
“’ OK gm ce-° 
0.00 0.00 e6>l 1,000 i684 1,000 1,28 18,418 
1.00 0,60 4 0.867 1677 0.880 1,08 18,126 
1.97 1,00 .749 0.064 1600 0.945 0.95 16,882 
2,98 1.50 s<655 0,486 Ilble 0.685 0,65 16,604 
5,87 2.00 «S87 0.666 1448 0.643 U.77 14,840 
7,46 4.00 140 0.751 1141 U,07 0.56 14,010 
15,68 6.00 1628 0,571 78 0.471 0,86 85s 
¥B.21l, 16.00 1145 0,401 S46 0.512 0,26 6.515 
(¥) (10) (11) (12) (13) (14) (15) (28) 
wits c# c Isp Isp/Isp° Cr h h/h® 
ft secm+ rt sec71 sec pi, 
0.00 5775 T8bz 245.8 1,000 1,36 458 1,000 
1,00 bys 8107 851.8 1.083 1.57 474 1,085 
1.97 5976 be54 265.7 1,049 1.38 466 1,022 
2.93 6012 6580 “58.4 1,060 1,36 458 1,000 
5.87 6044 BL65 e5u.8 1,066 1.36 455 0,983 
7.46 6095 8459 J6..7 1,076 1.39 409 0,893 
15,88 60835 8596 L60.7 1,069 1,58 3v6 0,712 
23,21 5886 81c8 w62,4 1,065 1,36 227 0,496 
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TABLS AL 
VARIATION Of PARPORMANCE PAR baTHR., / 12K WEIGHT PRE ChAT 
OF LIQUID slaUNin PORK HgO.-lHg(1)-N Hg SYSTEM 


ASSUSING nguLLIERIUN FLOW CONDITIONS 


(2) 


(o) 


(4) 


(8) 


(6) 


(7) 


(8) 












= — Ff, Be” » wae oa ¥ 
. x °*%K gm en™? 
0.00 0,00 ° 2852 1,000 1979 1,700 1.88 19,747 
4.08 0,25 “60% 0,983 1767 0.893 1,88 Si Jeea 
7.84 0,60 2677 «0,939 «891606 0,612 1.19 18,049 
14,54 4,00 2376 «0,883 «491542 «400.676 1,12 16,719 
20.3% 1,50 2099 0.7356 1140 0.576 41,07 15.696 
265,68 #00 1862 0,655 S77 0,494 1,08 14,801 
(¢) (10) (41) (Gk) 8) (14) (15) (16) 
WT cx e Isp Isp/Isp Cp oh h/n° 
£t sec7+ ft sec7? see mi. 

0.00 5773 8087 260.5 1,000 1.40 489 

4.58 5804 8055 wau.5b §60.96 1.88 479 

7,64 5699 TRE c46,0 0,982 1.38 460 0, 
14,54 D487 7648 257.5 0,946 1.39 418 0, 
20.32 5285 761 - 238.8 o.9ll 1.89 371 047 
25,38 5081 7048 216.9 0,874 1.389 $33 0,661 








A=-le 
TABLE XII 


VARLAT LUN Ue Part L ARCS PaRAR eT! Ao sith 


(1) 
WTp 


G00 
4,08 
7,84 
14,54 
20,52 
25,58 


(¥) 


0.00 
4.08 
7,84 
14.94 
20.d¢ 
xD. 58 


OF LI ylbin 


(<) 


BOLAS 


Q,00 
0,26 
J.50 
1.00 
1,50 
=.00 


(lu) 


ct 


AMUORIA FPUn HeOw-WRg(1)-Sel, SYS TEx 


Ao Ub ING Cu iu (a tel 


T COmPUsl TION FLOW COMLITIONS 


Bivit ruc CinT 


ft seem) ft secml 


(5) (4) (5) (8) 
;. 2a: >» *»%” 
O% 7 
1 1,000 1694 1.000 
e004 «=60. UHS—(<‘id CCL UTE 
“577 Q,985¥ 1555 4 =O. S17 
2576 0.855 1355 0,767 
2099 8 0.756 1138 )0— «0, 6 72 
1662 0.653 977 0.577 
(11) (lz) (15) 
c Isp Isp/Isp 
sec 
TBDE #45.8 1,000 
7316 =45.8 1.008 
7376 244.6 1,005 
7657 ESTs2 0,966 
7246 228.1 0, 956 
7048 218.9 0,488 


(7) (8) 
? Y 
go cu? 
1.28 18,416 
1.23 16,782 
1,18 16,601 
1,22 16,710 
1.07 15,686 
1.08 14,901 
(14) (15) (18) 
Cp bh h/k®° 
mi, 
1.86 460 1,900 
1.36 462 1.004 
1,38 458 0,986 
1.39 419 0,911 
1.39 #71 0,607 
1,39 333 0,724 


A~15 


TABLE XIII 


VARIATION OF PERFPORBANCES PaRAbhLTENS With #elGhi PiR CENT 


(1) 


WTS 


0,00 
0.90 


WI% 


0.00 
0,90 
1.79 


O,ol 
6.78 


15.46 


21,45 


ASSUBING BQUILIBHIUM FLOW CONDITIONS 


(2) 


MOLES 


(3) 
Te 


(11) 


C 


ft sec 
7870 
8078 
8107 
8160 
8193 


amy 
Bz 


8179 
7889 


(4) 
Ta/Te 


1.300 
0,997 
0,966 
0.929 
0,890 
Ue 704 
0,568 
0.400 


oO 


(©) 


(12) 


“1 


Isp 


See 


(6) 


0 
,/% 


i ,000 
0,912 
0,842 
0.782 
0.7351 
0.573 
0.401 
O.<76 


(15) 


Isp/Isp° 


244.4 


£43 ,5 


avl.& 


BOO 


~O4 .4 


256.1 


v4.0 
246.0 


1.0090 
1.0<0 
1.050 
1.057 
1,041 
1,048 
1,059 
1,00z 


(7) 
- 


gn cm 


1,28 
1,11 
0.98 
0.88 
0.79 
0,99 
0,40 
0.07 


(14) 


Cp 


O# LIQUID BHYDAUGaN FOR HPNA+Hy (1)-Nolg OYSToM 


oO 


(8) 
i 


1.367 
19.791 


(18,4138 


17,201 
16,153 
13,064 
9, 6Y8. 
7,019 


(15) (16) 


h/h? 


1,000 
1,Q17 
0.836 
0,976 
0,948 
0,849 
0,680 
0,461 


A-14 
TABL!. AIV 
VARIATION OF PERFORMANCE PAKARETEAS FITH “SIOKT PIR ChnT 
OF LIQUID HYIMHOGhN 70h RPNA-Hy(1)-NLR, SYSTEE 
ASSUMING CONSTANT COMPOSITION FLOW CORLITIOWS 





7 
(2) (2) (3) (4) (5) (c) (7) (8) 
7? O66M % Tete To tTe/Za F ¥ 

4 Oo, gs ca-> 
0.00 0.00 ebb? 1.000 1703 1.000 1.28 £#£1,006 
0,90 0.50 e945 U.6U7 l6v6 U.9¥6 1,1] 16.614 
1.79 1.00 2806 0,966 1618 0.951 0,86 16,820 
«,65 1.80 #746 0.928 1840 0,004 0,68 17,149 
$,51 #.00 265z 0,600 1655 6,554 0.79 16,126 
6.78 4,00 2B 0.784 1175 #«#+O,6868 0.56 18,068 

13,46 8.00 1660 0,568 616 400.478 0,440 » 683 

21.45 15,00 1184 # 0,400 645 0.319 «0,27 7,019 
(9) (10) (41) 9 (le) ~— (8) (a4) (25) (16) 
wip ce c Isp Isp/Isp® Cr h b/h® 

ft secm+ rt secm+ sec pi, 

0,00 50660 7656 7.2 1,000 1,85 437 1.000 
0.90 5811 7882 ca4.8 1,0u8 1.4 447 1,088 
1.72 5644 b0e< c4¥.1 1.046 1.57 447 1,025 
4,65 5698 as) L51.4 1,008 1.57 441 1.008 
3.51 8908 8155 &S5,3 1,065 1,368 450 9,984 
&.78 5958 8x04 266.0 1.075 1.38 S86 0,683 

12.46 5927 8179 294.0 1,066 1,38 315 0.716 

£1.44 »764 7868 445.0 1.030 1,36 212 6,485 


A-15 
TABLE XV 
VARIATION OF PERFORMANCE PARARSTARS WITH WEIGHT PER CENT 
OF LIQUID AMMONIA FOR RFNA-NHs(1)-NoHg SYSTEM 


ASSUMING BQUILIBRIURK FLOW CONDITIONS 


(4) 


(1) (2) (3) (5) (8) (7) (8) 

WI MOLES 1, #£Te/Te® Te Te/Te® i 
% °K gm on7° 

0.00 0.00 2957 1.000 2057 1,000 1.28 21.367 
3.70 0.25 2913 0.985 1823 0.685 1.24 +0,584 

7.13 0.50 2783 0.941 1648 0,809 1.20 16,491 
13.51 1.00 £479 0.838 1585 0.680 1.18 17.995 
18.72 1,50 2195 0.742 4177 0.5768 1,08 16.838 
35,49 2,00 1948 0.652 1008 6.495 1,04 15,923 

(9) (10) (11) (12) (13) (14) (15) (16) 
wre ci c Isp Isp/Isp? Cp 2b nh/b® 

ft sec ft sec sec mi, 

0,00 5660 7870s -k44,44 81,000 1,39 460 1,000 
3.70 5703 78863 8=6244.e «600858 )0=—ssa1,38 4460 1,000 
7.13 5887 7774 241.4 0.988 1.89 441 0,959 
15.31 6418 7508 «235.2 «20.954 ~=-1,.39 405 0,880 
18.72 5206 72c8 224.5 0,919 1,89 358 0.778 
23.48 5018 6858 216.1 0.884 1,389 314 0.683 


A~-16 


Tables AVI 


VAKIATIUN OF PUnPuUnBANCE PARAREYeuw ITH WelOwy Per CENT 


(1) 


0,00 
*.70 
7,15 

13,51 
“18.72 


oo. 4Y 


(3) 


0.00 
5.70 
7,15 
13.51 
18,72 
2.49 


OF Liybly AMMUNIA PUN KEN otha (1) -B Ry, oYoTie 


ADCUMING CONSTART COMPOSLTIUW FLOW CONDITIONS 


(10) 


ce 


vod 
L706 
5587 
0418 
0206 
9018 


ft sec f£t sec” 


(4) (5) (6) (7) (8) 
une” & %/%* Ff t 
O, gn caW 
1,000 1708 1,000 1,88 £1,006 
0.865 1668 0.6868 1.24 20,216 
O.e41 1576 0.827 1.20 10,414 
0.638 1363 0.8600 1,13 17.977 
O,?74e 1173 0.668 1,08 16,640 
0.859 1008 0.682 1,06 16.983 
(11) (12) (15) (14) (18) (16) 
¢ len Ipp/lep Sp & n/d” 
1 sec mi, 
7667 #381 1,000 1,35 437 1,000 
775. 40,1 1,008 1,86 441 1,000 
7690 258.6 1,003 1.88 454 0,#69 
7487 Rk32,5 0,976 1436 402 0,920 
7209 WBS.@ «= «HOS «a, 38 «354 «(0,610 
6958 Yl6,1 0.908 1.39 314 0,719 


a-lL? 





TABLE AVIT 


VARIATION OF B(T.)) SIYH WelGHT PAK Conv OF LIQUID EYDROGLN- 
FOR EQUILIBRIUR FLO“ CONDITIONS 


Propeliunt systems; 
A, NO,-H,(1)-NLUH, 
B. Hs0y-H, (1)-N Hy 
Ce RPNA-H, (1)-HHy 


Wt.% Ho(1l) A B C 
0.00 0.v00 3.000 0,090 
3,00 0.108 0.110 0,124 
6,00 0,236 04240 0,26 
8,00 0.520 0.414 0.329 
10.00 0,385 0,574 0,388 
12,90 0,440 0,428 0.436 
15,00 0,506 0,496 0,504 
20,00 “ 0.580 0.587 
Averaze Error 0.004 0,006 0,006 
Max, Error 0,008 0.008 0.010 











Errors are with respect to average value obtained 


by averaging B (155) of systems A, B, and C, 


AWji8 
TAELi. AVITI 


VAHIATION OF "(Igp) WITH WELUNT Per CANT Cy LIVUIL WY UROGLN 
FOR CONSTANT COMPOSITION *LOn CUALLTTIONS 


Propélliant Systems: 
A. W, Og-Hy (1) -N,Hy, 


C. RPNA=Hy(1)-N,H, 
D 


. 05 (1)-Hy (2) -H ny 


Average 

We.g H.(1) 5 B ¢ of A, Bas D 
0.00 0,000 0.0% 0,000 0,000 0.000 
5.00 0.114 0.14 0.144 0.161 0,14 
6.00 0.260 0.<06 Q0.«8L 0.<76 0.<50 
8.00 J, 565 UO, dao 0, 508 0,056 0 oa 
10,00 0.456 0.400 0.416 0.417 0, 582 
1z.00 0.440 0 .abz 0.464 0.469 0,452 
1,00 QO, 056 0, o&e 0,550 0,556 0,638 
20.00 - 0,606 0.614 0.610 0,630 
Average errord.0ls 0.011 U.005 0.016 
Max, trror 0,0<1 0,017 0,013 0,0u6 


trrors are #ith respect to Average value obteined 


by averaging E(Isp) of systems A, Band C, 


(a) D&ta was computed from values listed in hef, 3, 


A-19 





TABLE XIX 
VARIATION OF E(py WITh WEIGHT Pe CENT OF LIQUID HYDROGEN 
FORK EQUILIBRIUm FLOW CONDITIONS 


Propellant Systems; 
A. Ny04-Hz (1)-NzHa 
B. H20g-H2(1)-N2Bg 
C, RFNA-He(1)-livlig 


Average Error 


Max, Error 


0.006 


0.007 


0.006 
0,016 


Wt. Ho(1) A B c 
0.00 0,000 0.000 0,000 
3,00 0,040 O.0e# 0,052 
6.00 0.080 0,064 0,074 
8.00 0.004 0,088 0,080 

10.00 0.100 0,094 0,089 
12,00 0.105 0,095 0,096 
15.00 0.108 0,096 0.100 
20,00 ~ 0.080 0.068 


0,006 
0.014 








Errors are with respect to average value ob- 


tained by averaging (hh) of systems A, B and C, 


A-O 


VANLATION OF E(,) WITH WEIGHT PLA CeNT OF LIQUID WYLnOGeH 
FOR CONSTANT COMPOSITION FLOW COW TIOAG 





Propellant Systems; 

A. Nz0g-Hy (1) -Woky 
B. HpO,-H, (1)-Wpky 
C. RPMA=|Ey (1) =H, 











by everaging E(h) ror systesis A, B anu C, 


‘Wt.e He (1) A B c a o 
0,00 0.000 0,000 0,000 0.000 
3,00 0.113 0,080 0,086 0,083 
6,00 0.164 0.150 0.116 0,137 
8,00 0.184 0.149 0,128 0.184 

10.00 0.188 0,158 0,156 0.169 
0.190 0,148 0.143 0.160 
0.180 0.144 0.140 0,155 
a 0,124 0.106 0.115 
Average Error 0,027 0,009 0,017 
(Max, Brror 0,050 0.018 0,026 


Errors ure with respect to overage Vvilue obtained 


A-zl 
TABL:| AAT 


VARIATION OF E(1,)) WITH WEIGHT PEK CENT OF LIQUID AMKONI 
FO EQUILIBRIUk FLOW CONDITIONS 


Propellant Systeus; 
D. W,0g-NHy(i)-N2Hy 
B. HL O.=NHy(1)-NUHy 
F, RYNA-NHz(1)-NUHg 


Errors are with respect to average Value obtainea by 


Wt. NHz(1) D E F 
0,00 0,000 0,000 0,000 
3,00 0,022 0.010 0,010 
6.00 0,03x 0,0x8 0,034 
8.00 0.044 0.044 0,056 

10.00 0.060 0.064 0.078 
12,00 0.078 0,088 0.100 
15,00 0,106 0.120 0,136 
20.00 < 0.172 0,190 
Average Error 0,007 0.00% 0,008 
Max, Error 0,015 0,009 0,015 


sweraging B(z),) for systems LD, E ana F, 










TABLA KXII 


VARIATION OF (Igy) #ITM WEIGHT PER CENT OF LIQUIL aenowta 
a FOs UvWeTaNt COMPOSITION FLOW CONUITIONS 





ropelient Systems: 

D, W.0q-WH), (1) -Np Hy 
Be A. 0, -WH, (1) 8 Hg 
FP. APMA-Mg (1)-NyHy 


’ 




















yt. Whig (1) D : P 
0,00 2.000 0.000 0.900 
8,00 0,020 0,.0%0 0.u18 
6.00 0,044 0.044 0.050 
6. 90 0,063 0.Uu64 U.076 
“0 co 0.004 0,086 0.100 
14.00 0.107 0.110 0.123 
— 45,00 0.141 0.144 0,158 
80.00 - 0.186 0.210 
rage Error 0,004 0,003 0,007 
Mex, krror 0.006 0,007 O10 





wn, 





Errora are with respect to average velue obtained 
by avervging L(Igp) for systems b, E, and F, 





A-~<d 


TABLE XAIIT 






VARIATION Of Bp) WITH WhIGKT PRR Cert OF LIQUID AKMONIA 
FOR EQUILIBRIUM FLOW CONDITIONS 


Propellant systens: 
D. WoOg-NH,(1)-NUHy | 
BE. HoOg-NHs (1) -Nikg s 
F.  RFNA-NB;(1)-NoHg | 


e a 





averaging E(n) for systems D, E, ana F, 


Wt.% NHs (1) D E i D, E& F 
0.00 0.000 0.000 0.000 
3.00 0.012 0.000 0,010 
6.00 0.020 0.000 0.016 
8.00 0.020 0,003 0.0288 

10,00 0,0Z1 0.008 0.052 

12,00 0,024 0.016 0.040 

0,028 0.02c8 0.042 

20.00 - 0.028 0.08% 

Average Error 0,004 0,009 O.007 
Max, Error 0.008 0.012 0.0135 


Errors are with respect to average value obtained by 








A-<4 
TABLE KALV 


VAKIATION OF L(y) WITH WEIGHT PER CENT OF LIQUID AMMONIA 
FOR CONSTART COMPOSITION FLOW CONDITIOSS 





Propellant Systems; 

D. NzOg=NHg (1) =iyHy 
B, HyOp=Nhig (1) -tiy 
F, RFNA=NHy (1)-yH,g 





D - r estes F 
0.00 ¥.000 0.000 0,000 0.000 
©,00 0,026 0,016 0.018 0,021 
6,00 0.060 0,042 0,042 0.048 
8,00 0.076 0.0668 0,060 0,065 
10,00 0,088 0.070 0,072 0,077 
12,00 0.10% 0,080 0.080 0,088 
15,00 0.106 0,060 0,064 0,080 
20,00 - 0.070 0.064 0,067 
Average Error 0,012 U.007 U. 008 
Max, &rror 0.016 0.ULO 9,008 








Errors are with respect to «veruge value ob- 


tained by averaging f(n) for systems DP, B ang F, 
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